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Abstract. The multigraded Poincare-Betti series t) of a monomial ring 

k[x]/(M) on a finite number of monomial generators has the form IIic g„(l + 
X it)/°R fc(»;t), where bji fc(x;t) is a polynomial depending only on the mono- 
mial set M and the characteristic of the field k. I present a computer program 
designed to calculate the polynomial £>u,fc f° r a given field characteristic and 
a given set of monomial generators. 



1. Introduction 

Let Q — k[»] = k\x\, . . . , x r \ be the polynomial ring over a field k with r variables. 
The ring has a natural INT-grading 1 by setting deg(xi) = for the canonical basis 

vectors et of N r . Setting |(oi, . . . , a r )\ = ai H h a r , we can derive an N-grading 

of Q from this N r -grading. We write deg(m) for the N r -degree of a monomial m, 
and given (ai, . . . , a r ) = a S N r , we write a;" = a;" 1 ^^ 2 • ■ ■ i° r 

The N r -grading and its inherent N-grading both are inherited from Q to the ring 
R = Q/I where J is a monomial ideal, i.e. an ideal generated by monomials in Q. 
We call R a monomial ring. 

The N r -grading inherits, via minimal resolutions respecting to the grading, to an 
[NT-grading on Tor^(fc,A:). Thus, we can define the multigraded Poincare-Betti 
series Pj^fa t) of an INT-graded ring R over k: 



J2 J2 dim fc Tor f ik,k) a x a t l 



From Pj^(x;t), the simple Poincare-Betti scries can be calculated as Pj^(t) = 
Pf(l,...,l;i). 
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Since Backclin , 1982| , it is known that the Poincare-Betti series of a monomial ring 
R holds the form P*(t) = for some polynomial &R,fe(i). 

Moreover, Alexander Berglund proved in Bcrghm d, 2004| that deg < 2n for 

monomial rings with n monomial generators. It follows that there are finitely many 
Poincare-Betti serie occurring for a fixed n at all. 

We may define a partially ordered graph, or po-graph, to be a graph with a par- 
tial order on the vertices. Two po-graphs are said to be isomorphic if there is a 
simultaneous isomorphism of the graph and the partial order. 

Let Q = k[x] and Q' = k[x'] be polynomial rings for two finite variable sets x = 
{x\, . . . , x r } and v! = {x%, . . . , x r >}, and let / and I' be monomial ideals in Q and 
Q' respectively, with M and M' the sets of generators for each ideal. For some set 
S of monomials, we denote by L$ the set of all least common multiples of subsets 
of S. Ls can be equipped with the structure of a po-graph, ordering monomials by 
divisibility and adding an edge between two elements when they have a non-trivial 
common factor. 

Luchezar Avramov shows in |Avramov7 "2002| that if two rings R = Q/I and R = 
Q'/I', with / and I' generated by the monomial sets M and M' respectively are 
such that L M = L M ,, then b Ryk (t) = &R',fc(t). From this follows that for a fixed field 
and a fixed number of monomial generators, only finitely many different Poincare- 
Betti series can occur. Avramov further proves that the limitation to a fixed field 
is superfluous. 

In |Berglund, 2004| , Alexander Berglund, proving the conjecture by Avramov that 
deg(6i{,fc(i)) < 2n whenever the monomial ideal generating the ring R has n genera- 
tor, constructs combinatorially a minimal model for R and gives a characterisation 
of the Poincare-Betti series denominator polynomial in terms of the homology of 
associated simplicial complexes. Avramov's earlier observation that only finitely 
many different Poincare-Betti series exist for a fixed number of generators for the 
monomial ideal I follows as an immediate consequence of Berglund's construction. 

I have in the course of my M.Sc. thesis work |Johansson, 2004| continued Berglund's 
work by implementing his formula in a program, written in C++, capable of calculat- 
ing simplicial homology over fields of arbritrary characteristic, as well as explicitly 
calculating bn,k(t) for arbritrary characteristic of the coefficient field k. 



1.1. Simplicial Complexes. Since we will work a lot with simplicial homology, I 
will take a few moments to review definitions and terminology. A simplicial complex 
on a set V is a set A of subsets of V such that if G € A and F C G, then F € A. 

V is called the vertex set of A. All simplicial complexes I shall refer to will have 

V = (J A unless otherwise stated. The i-faces or i-simplices of A are precisely the 
elements in A of cardinality i + 
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To a simplicial complex A we can associate an augmented chain complex C(A) 
with Cj(A) the free abelian group on the i-faces of A, where we consider to be 
the unique — 1-simplex. We equip C(A) with the standard differential of degree 
-1. Thus 

Hi(C(A)) = Hi(A) 

As usual, for an abelian group G, we set C(A; G) = C(A) ®j_ G and ifj(A; G) = 
Hi(C(A',G)). 

For a simplicial complex A we define the Alexander dual 

A v = {F C V | V \ F <£ A} 

For simplicial complexes A and A' with disjoint vertex sets we define the join 
A * A' = {F U F' | F e A, F' e A'} 

and the dual join 

A • A' = (A v * A' v ) v 

and note that then 

A ■ A' = {F e V U V | F n V G A or F n 7' G A'} 

In |Bruns and Herzog, 1993| Lemma 5.5.3], it's shown that 

ff i (A;A;)^# n _3- i (A v ;fc) 
for a complex A on n vertices. 

For a graded vector space H = ^ iei4 Hi, we will write H(t) for the generating 
function J^ieM ^ m k Hit 1 of H. We then can find that the join of complexes as well 
as the Alexander dual gives rise to rather easily handled equalities on the level of 
generating functions for their respective homologies, again with A a complex on n 
vertices, and A' some other complex. 

(1.1) f\ff(A v ; = t 3 H(A; k){t) 

(1.2) H(A * A'; k){t) = tH(A; k)(t) ■ H(A'; k)(t) 

where the t factor in the latter equation comes from the fact that a simplex with d 
elements is considered to have dimension d—1, 

2. Berglund's work 

In |Berglund, 2004| , Berglund treats the theoretical aspects of computation of the 
Poincare-Betti series denominator for monomial rings. I will not repeat all of his 
proof here, but rather reference his work to establish the vocabulary and touch the 
results I will need for my own work. We can find Q = fc[x] is [NT-graded by assigning 
to a monomial x^ 1 . . . x^ r the element a = {a\, . . . , a r ). We write x a for a;" 1 . . . x% r . 
The monomial x a is said to be squarefree if a G {0, l} r . 
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By a construction by Weyman and Froberg W eyman, 1977[ |Froberg, 1982| , it is 
enough to treat squarefree monomial sets, since an easy procedure can be used to 
go from a monomial ring to a squarefree monomial ring with the same homological 
properties. So we can assume that R = Q/I is a squarefree monomial ring with / 
generated by the monomial antichain M of cardinality n. 

For a set S of monomials, set ms — lcm(m : m e S). In particular = 1. Thus 
Lm = {ms | S C M}. For a specific monomial m and a monomial set M set 
M m = {m! e M | m'\m}. 

Now, for / = (M) a monomial ideal in Q generated minimally by the antichain M, 
Berglund introduces the complex 

Am = {S C M | ms ^ wim or S disconnected} 

where connectivity is for S as a subgraph of Lm- 

Using multigraded ring-deviations 

e iia = dim fc -Hi_ 3 (A AfxC , ; k) 
Berglund gives the squarefree part of the multigraded Poincare-Betti polynomial 

(2.1) b R , k {x,t)= Yl {l-x a p a {t)) (mod (xl,...,x 2 r )) 

with p a (t) = i 3 i/(Aj\,f xCt ; k)(t). Backelin demonstrated already in |Backelin, 1982| 
that the denominator polynomial will be squarefree whenever the monomials gen- 
erating the ideal all are. 

Berglund then goes on to find several more theoretically pliable forms of this par- 
ticular formula; expanding the product and taking the irrelevance of non-squarefree 
terms into account, he arrives at the form 

(2.2) b R , k (x; t)=l + J2 m s (-t)^ s}+2 H(A s ; k)(t) 

s 

where the sum is taken over all non-empty saturated subsets of M and c(S) counts 
the number of graph components of S as a subgraph to Lm- We define the satu- 
ration of a subset S C M as the set of all monomials in M that divide the least 
common multiple of some connected component of S as subgraph to Lm- A set is 
saturated if it is equal to its saturation. 



3. The resulting application - poincare 

My own achievement is that I have constructed a computer program to calculate 
simplical complex homology and mainly to calculate the denominator polynomials 
of Poincare-Betti series using Berglunds methods. I will devote this section to a 
discussion of the program, which can be fetched in its latest version under the MIT 
software license from http://www.math.su.se/~mik/poincare/ 
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The form deemed most promising for implementation as I started was the form 
given in l|2.1|l - mainly since the formulation in terms of saturated subsets had 
at that time not yet matured. Thus, I have implemented specific C++ classes for 
calculating in the ring Q/(xf, x%, . . . , x 2 .) and let the final product forming the 
polynomial take place in that particular ring. The only part forming any kind of 
complexity for the straightforward implementation is that of forming the complex 
Ajh i0 and calculating its homology over the specified characteristic. 

The construction of Am xC « is done with a modified kind of breadth-first search: 
monomials are stored in a queue along with an index keeping track of which of the 
monomials covering the particular monomial that have already been tried. Thus, 
for each monomial in the queue, all later covering monomials are tried one after 
the other, and upon compliance with the two conditions - that the least common 
multiple of all generating monomials dividing the candidate is equal to x a and that 
those generating monomials are connected as a graph - the monomial is added to 
the queue carrying a testing index one higher than the index that produced it. This 
algorithm does yield a speed increase compared to the earlier algorithm that simply 
tested all monomials for both conditions; but still is not optimal by far. 

Once the simplicial complex as such has been constructed, the calculation of its 
homology commences. This is calculated degree by degree, constructing a matrix 
with entries in {0, ±1} and fetching its rank from the matrix routines in the Pari 
library [Pari/GP, 2004] . 

This is wrapped in a text-mode user interface, using the GNU Readline library 
to facilitate command history and command editing. The user interface reads in 
space-separated lists of monomials as input to the add simplex and add monomial 
commands. A monomial, to the program, is a ^-separated list of strings of charac- 
ters, where each separated string is taken to be the name of a variable. The variable 
names must avoid +-*/" , . and whitespace, but can use any other characters. Any 
string occuring in such a position will be interpreted as a variable and added to an 
internal dynamic variable pool. 

The user interface wraps, among other things, around the Weyman-Froberg method 
for conversion to a squarefree monomial ring. The conversion is done transparently, 
using several internal variables that are easily converted back to the original vari- 
ables before printing the answer. The output rendered by the program is written in 
such a way that other computer algebra systems should have an easy time handling 
it. 

As an example on how the program works, I give in example l3.H a session, calculat- 
ing first the simplicial homology of the projective plane over Q as well as over Z2 and 
then calculating the Poincare-Betti denominator polynomial of y, z]/ (x 2 , xy, yz). 

Among the things we may observe in example 13. II is the basic set of commands - 
add simplex and add monomial to build simplicial complexes or monomial ideals, 
homology and denominator (or an abbreviation denom thereof) to calculate sim- 
plicial homology and the Poincare-Betti denominator respectively, as well as the 
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Welcome to the Poincare calculator. You can use this program to 
calculate simplicial homology over prime fields and to calculate 
the denominator polynomial of the Poincare-Betti series of monomial 
rings . 

(c) 2004 Mikael Johansson 

This program is released under the MIT License 

> add simplex a*b*e a*b*f a*c*d a*c*f a*d*e 

> add simplex b*c*d b*c*e b*d*f c*e*f d*e*f 

> homology 

Calculating homology ranks . . . 

***** Hilbert series of simplicial homology ***** 


> char 2 

New characteristic: 2 

> homology 

Calculating homology ranks . . . 

***** Hilbert series of simplicial homology ***** 
ZZ + ZZ~2 

> add monomial x~2 x*y y*z 

> char 

New characteristic: 

> denom 

1 - x~2*ZZ~2 - x*y*ZZ~2 - y*z*ZZ~2 - x~2*y*ZZ~3 - x*y*z*ZZ~3 

> set multigrade false 

> denom 

1 - 3*ZZ~2 - 2*ZZ~3 

> quit 

Thanks for visiting. 

Example 3.1. Calculation with poincare 



command char, which changes the field characteristic over which all homology 
calculations take place and the command sequence set multigrade false, which 
sets a flag that causes the program to change the way it prints the polynomials out- 
put by the denominator command, so that instead of the polynomial &R,fc(x, i) the 
program prints the polynomial . . . , 1, t). The program ends upon receiving 

quit. 

In addition to these, there are the commands clear: clearing the stored simplicial 
complex and monomial ideal, but not changing the characteristics used and var: 
which changes the implicit homology variable, which in my review in this paper has 
been called t, and which by default in poincare is called ZZ. Should you wish to 
use ZZ as a ring variable in your calculations, a change of homology variable will 
be necessary. For this, the command var is provided, with which you can change 
the string that poincare uses for the homology variable. 
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4. Questions and future directions 

There are several things that I want to improve upon on the system herein pre- 
sented, and also several questions that can be posed. 

There are numerous complexity issues associated to the program in its current 
form. Mainly, these issues are related to the size of the resulting po-graphs for 
larger sets of monomials. As an example, the initial ideal of a Grobner basis with 
revlex ordering of the homogenized cyclic 6-root ideal, i.e. the ideal generated by 

xi + x 2 + x 3 + x A + x 5 + xq, 

xix 2 + x 2 x 3 + x 3 x 4 + x 4 x 5 + x 5 x 6 + xix 6 , 

X\X 2 X 3 + X2X3X4 -I- X3X4X5 + X4X 5 X 6 + XlX 5 X 6 + XlX 2 X 6 , 

XlX 2 X 3 Xi + X2X3X4X5 + X 3 X4X 5 X 6 + XlX 4 X 5 X 6 + XlX 2 X 5 X 6 + XlX 2 X 3 X 6 , 
X\X 2 X 3 X 4 X^ + X 2 X 3 X\X$X§ + XiX 3 X4X$X6+ 
XiX 2 X4X 5 X & + X 1 X 2 X 3 X 5 X 6 + XiX 2 X 3 X4X 6} 

xix 2 x 3 X4X 5 x 6 - y 6 

invariably becomes larger than Pari's working memory. This initial ideal has 100 
monomial generators and produces a po-graph with 11443 elements. The calcula- 
tions normally halt after between 200 and 400 lattice point calculations. 

The problem I have observed with for instance this example is that the calculation 
of homology of large simplicial complexes is memory wise unfeasible. The most 
visible problem is when the homology calculations turn out to be too hard, since 
this results in a crash in the Pari library; whereas too hard construction of the 
simplicial complexes merely result in slow running of the program. 

Question 4.1. Given the rather special structure of the matrices that are used 
to calculate field homology of a simplicial complex, can anything be said about the 
sizes of elements of a matrix as some sort of reduction algorithm is used to deduce 
its rank? Can some variant of row- reduction be found such that the matrices that 
produce "large" entries (for instance, larger than a standard 32-bit word) have some 
easily recognizable feature ? 

If such a method could be found, then it would be possible to single out the specific 
homology matrices that would actually need treatment with some sort of bignum 
library, and produce fast 32-bit (sparse) matrix arithmetic to deal with all other 
instances. Such a separate treatment would also obliterate the need to rely on 
Pari's internal memory allocation structures and would give the user more control 
over acceptable memory consumption for the calculations at hand. 

An obvious further point of attack is the generation algorithm for the Berglund 
complexes. This would be vastly improved if a search algorithm would be con- 
structed that minimizes the number of multiple checks done on each candidate 
monomial; since the graph connectivity checks are not, in the context, very fast. 
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Finally, an interesting direction to take would be to look at the APIs for larger 
computer algebra systems and try to adapt the code here written to work as a 
pluggable module to those systems; for instance providing an interface to calcu- 
late Poincare-Betti series of monomial rings directly from Singular or Magma or 
Macaulay 2. 
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